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ABSTRACT 

Aims. We studied the global distribution and kinematics of the extra-planar neutral gas in the Milky Way. 

Methods. We built 3D models for a series of Galactic Hi layers, projected them for an inside view, and compared them with the Leiden- 
Argentina-Bonn 21 -cm observations. 

Results. We show that the Milky Way disk is surrounded by an extended halo of neutral gas with a vertical scale-height of 1.6^^ 4 kpc and an 
Hi mass of 3.2^^^ x 10^ M©, which is ~ 5-10% of the total Galactic H i. This H i halo rotates more slowly than the disk with a vertical velocity 
gradient of -15 + 4 km s~^ kpc~^ We found evidence for a global infall motion in the halo, both vertical (20^^ km s~^) and radial (30^^ km s~^). 
Conclusions. The Milky Way H i extra-planar layer shows properties similar to the halos of external galaxies and it is compatible with being 
predominantly produced by supernova explosions in the disk. It is most likely composed of distinct gas complexes with masses of ~ lO'^'^M© 
of which the Intermediate Velocity Clouds are the local manifestations. The classical High Velocity Clouds appear to be a separate population. 
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1. Introduction 

In the recent years it has become clear that neutral hydrogen 
in nearby galaxies is not confined to a thin disk, but it extends 
for several kiloparsecs into the halo region. Several processes 
may contribute to create this halo (or extra-planar) gas: galac- 
tic fountains and extra-galactic accretion are the main candi- 
dates. In the first case the gaseous halo is built up by the super- 
nova feedback. In the classical picture f Shapiro & Field|1976t 
|Bregman,19 80 ) supemovae explosions heat the gas in the disk 
and eject it up to a distance of few kiloparsecs, the gas travels 
through the halo region, it cools and eventually falls back to the 
plane. In the latter scheme the halo originates from gas accre- 
tion out of the intergalactic environment ( |Qort 197Q[ Kaufmann 
|et al.||2Q06] ). The chemical and kinematical properties of the 
halo gas expected in this second case can diff'er significantly 
from the properties of the gas in the Galactic disk. 

Observations of external galaxies suggest that both pro- 
cesses contribute. A large amount of halo gas has probably an 
internal origin. High velocity outflows of neutral and ionized 
gas from the disk are observed to be correlated with star form- 
ing regions ^Kamphuis et al.|19"9T Fraternali et al.|20()4 l. The 
extra-planar neutral hydrogen has a kinematics similar to that 
of the disk H i ( [Fraternali et al.||2002[) but with some pec uliar 
properties: a) it rotates more slowly ( Qosterloo et al.|20Q7| ); b) 
it shows vertical motions ( [Boomsma et al. 2005 1 ); c) it has a 
global inflow ( [Fraternali et al.|200T] ). Extra-planar difl'use ion- 
ized gas (eDIG) is present in all galaxies with sufficiently high 
star formation rates ( [Rossa & Dettmar[2003| ) and its kinematics 



seems to follow closely that of the neutral component ( Heald 
et al.|2006[[Kamphuis et al.|2007j ). 



A fraction of the halo gas may also have an external ori- 
gin: some galaxies (including the Milky Way) are surrounded 
by clouds or filaments of H i with highly anomalous kinemat- 
ics, as in the cases of M 3 1 ( [Thilker etal .|2004|) and NGC 2403 
( [Fraternali et al.|20Q2] ), that may be the signature of gas accre- 
tion and minor mergers. This extra-galactic gas accretion is of 
vital importance to sustain the star formation in galaxies like 
the Milky Way ( jSancisi et al.|2008| ). 



In the Milky Way the Intermediate and High Velocity 
Clouds (IVCs and HVCs, [W^er & van Woerden|1997| ) con- 
stitute the main evidence for the presence of extra-planar gas. 
The HVCs are clouds of neutral hydrogen with a line-of- sight 
velocity significantly anomalous with respect to the Galactic 
rotation. They have been known for half a century ( [Miinch 
& Zirin|ri961| ) but only in the last decade observations have 



unequivocally confirmed that they belong to the halo region 
( Wak ker|[2Q5^ jWakker et aL][200gl ) at distances from the Sun 
< 10 kpc. Also, they are observed to have low metallicity 
(-0.1 Zq), which suggests an external origins for most of these 
objects (see for instance [Tripp et al.|2003 [ for complex C). In 
contrast, the IVCs (H i clouds with a lower anomalous veloci- 
ties) seem to constitute a diff'erent population, since their dis- 
tances are lower (~ 1 kpc), their metallicity is higher (~ 1 Z©) 
and their kinematics follow more closely the Galactic rotation. 
Therefore these clouds have been regarded as possible galactic 
fountain objects ( Wakker|2001 ). 
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Anomalous H i clouds are not the only observable compo- 
nent of the Galactic halo. A large amount of diffuse Hi emis- 
sion has been found by ILockmanl (| 1984|) in the lower halo re- 
gions of the inner Galaxy; Levine et al. ( 2QQ8| ) analysed this 
gas within 0.1 kpc of the mid-plane deriving a vertical lag 
in the rotation velocity of -22kms"^ kpc~^ Signs of a local 
halo building from recent stellar feedback are visible in the gi- 
gantic Ophiucus superbubble, this is a huge structure of neu- 
tral and ionized hydrogen reaching heights of ~ 3.4 kpc above 
the disk ( [Pidopryhora et al.||2Q07| ). The presence of an extra- 
planar DIG layer in the Galaxy has been known since long 
time (" Reynolds et al. 1973| ). Measurements of O vi absorp- 
tion lines in the Galactic halo obtained with the Far Ultraviolet 
Spectroscopic Explorer ( Sembach et al.|[20Q6 ) confirmed the 
presence of a 'hot Galactic halo' extended over tens of kilopar- 
secs, already hypothesized by |Spitzer| ( [1956| ) as a medium to 
pressure-confine the HVCs. 

In the literature there is little attempt to describe the Hi 



halo of the Milky Way as a global structure. Kalberla & 
|Dedes| ( |2QQ8l hereafter KD08) described the whole Hi distri- 
bution in the Galaxy via a spatial reconstruction of the Leiden- 
Argentine-Bonn (LAB, Kalberla et al. 2 005 ) H i Survey dataset. 
The de-projection of the data relies sensitively on the rota- 
tion curve assumed, which has been derived by modelling the 
Galactic dark matter halo ( [Kalberla et al.|2007] ) and under the 
assumption that gravity dominates the dynamics of the gas ev- 
erywhere. KD08 found that 10% of the Galactic Hi is located 
outside the disk and it is highly turbulent. 

In this paper we model the global H i component that sur- 
rounds the disk of our Galaxy. We make no assumption about 
the dynamical state of this gas but simply derive its global prop- 
erties directly from the data. Section |2] describes how we mod- 
elled the halo emission and what is the eff'ect of the various 
kinematical parameters. In Section [3] we derive the vertical ro- 
tational gradient in the inner halo and in Section]?] we estimate 
the kinematical parameter of the Galactic halo. A discussion 
of the physical interpretation of these parameters follows in 
Section |5l 



2. Modelling the Hi halo 

We assumed that the Galactic H i emission is due to two distinct 
components: the disk and the halo. The physical and kinemat- 
ical properties of the first are assumed to be known (Section 
2.1 ). In contrast, the halo properties are not known and its dis- 
tribution and kinematics have to be derived by fitting a number 
of free parameters. We modelled the Galactic Hi emission by 
building pseudo-datacubes with the same resolution and total 
flux as the LAB dataset, and we compared them with the latter. 

In the following we will use two coordinate systems. The 
first is the Galactic coordinate system, where / is the Galactic 
longitude, b the Galactic latitude, and vlos is the velocity along 
the line of sight. The second is a cylindrical system centred 
at the Galactic Centre, where R is galactocentric radius, the 
azimuthal angle (set to 0° at / = 0°) and z is the height from the 
mid-plane. 



2.1. Removing the H i disk emission 

In order to better visualize the extra-planar H i emission in the 
Milky Way we first removed the emission associated with the 
disk applying the standard technique that makes use of the de- 
viation velocity ( Wakker fl991[ ). For a fixed (/, b) pair, the de- 
viation velocity vdev is the diff'erence between the observed 
Vlos and the largest velocity predicted by a disk model. It de- 
pends on the rotation curve and on the disk geometry adopted. 
We assumed circular orbits and a flat rotation curve, indepen- 
dent of z, with rotation velocity v^(R,z) = Vq = 220kms~^ 
For the disk geometry we used the flared disk described in 
KD08, with the scale-height given by h,(R) = hQe^^'^^^'^' , 
where /z© = 0.15 kpc and Rf = 9.8 kpc. The hs(R) holds for 
R > 5 kpc and it is constant for inner radii. 

The line-of- sight velocity of a cloud at a distance d from 
the Sun in a generic direction (/, b) is given by: 



(R \ 
v^(R, z)^ - VqJ sin(/) cos(b) 



(1) 



where 7?© = 8.5 kpc is the distance of the Sun from the Galactic 
Centre. The maximum and the minimum values of eq. ([T]) (the 
'envelope' of vlqs) have to be determined by varying d in 
the range [0, Jmax], where Jmax is <^(^max,^max)- We assumed 
^max = 35 kpc (as in KD08) and Zmax is chosen to be 7 x hs(R), 
i.e. about 1 kpc in the central parts. With these assumptions the 
envelope is large, ensuring that the disk emission is completely 
removed. Given this envelope, vdev can be easily estimated for 
any (/, /?); in agreement with the literature ( |Wakker|| 2004 ), we 
cut out data with |vdevI < 35 km s"^ We also removed the Hi 
emission of the Outer Arm region (a distant and bright spiral 
arm that belongs to the warp; |Habing|]T966| ), the Magellanic 
Clouds, the Magellanic Stream and external galaxies. We stress 
that the geometrical parameters used to describe the shape of 
the disk weakly afl'ect the deviation velocity calculation. For 
example, assuming a disk with no flare and 2 kpc of thickness 
gives an envelope of vlqs similar to one above. 

Figure [T] shows the resulting velocity-integrated H i emis- 
sion obtained after the disk exclusion. All the H i emission in 
Fig. [T] is not compatible with a thin disk in difl'erential rota- 
tion. Some HVCs and IVCs can be recognized in this map, 
and they appear embedded in a more difl'use extended medium. 
This medium is visible in almost any direction of the sky and it 
is the main subject of our investigation. 

2.2. The H i haio distribution and kinematics 

We modelled the H i halo of the Milky Way assuming a density 



distribution p{R, z) analogous to that used by |Qosterloo et al.| 
P007| ) for the nearby galaxy NGC 891, described by the for- 
mula: 



p(^,.)=Po|l + ^) exp(-^)sech^(^ 



(2) 



where R^ = 1.61 kpc and y = 4.87, as recently estimated by 
Marinacci et al . | ( ]2"0 1 Ob| ) . For a given kinematics, the only pa- 
rameter that critically influences the resulting pseudo-datacube 
is the scale-height, while the scale-radius Rg has a minor efl'ect. 
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Galactic latitude 



Fig. 1. All-sky map of the Hi halo of the Milky Way: velocity integrated H i map obtained from the LAB Survey after removing 
the disk emission, the Magellanic Clouds and Stream and the Outer Arm. The sharp edges close to the mid-plane are artifacts of 
the disk removal. The dataset has been smoothed to 1° resolution. Contour levels are: 0.5, 10, 50 and 120 x 10^^ cm"^ in column 
density. 



Therefore we set hg as a free parameter, and we fixed Rg and y 
to the above values. Also, we set the central density po so that 
the total flux of the model equals the flux of the LAB data, once 
the disk emission has been excluded. 

We modelled the kinematics of the halo using three param- 
eters: 

1. the vertical gradient of rotation velocity dV(p/d\z\, so that 
v^fe) = V0 -h \z\ X dvcp/d\z\ for every R; 

2. the velocity along the z-direction v^, assumed positive if the 
gas is escaping; 

3. the velocity along the R-direction vr, assumed positive if 
the gas is moving to larger radii. 

At first, we assumed that these parameters do not depend on 
R and z. The resulting velocity field v = (vR,V(p,Vz) has to be 
projected along the line of sight to obtain vlos • 



VLos(hb,d) = vr 



R^ + d^ -Rl- 
2Rd 

(R, 



-h Vz sin(/?) + 



{^^0 - v© j sin(/) cos(/?) (3) 



where d is the heliocentric distance. Given p{R, z) and v(7?, z), 
the H I column density A^hi(^, v) or the brightness temperature 
T^{1, v) can be easily derived. We also included an isotropic 
velocity dispersion of 20kms~^ for the Hi halo, spreading the 
computed components Nmihb^v) over the neighboring chan- 
nels. 

2.3. Toy models 

To better understand the eff'ects of the diff'erent kinematic pa- 
rameters on the pseudo-datacubes we first built some 'toy' 
models with simple kinematics. For simplicity we set the Hi 
halo scale-height to 2kpc and its mass to 5 x 10^ M©. We re- 
moved the disk emission as discussed in 2. 1 Figure|2]shows the 



longitude- velocity slice (l-v) at the galactic latitude b = 30° 
for the LAB data and for six diff'erent models. The data panel 
shows that after the disk cut in vdev (white region at low ve- 
locities) a large fraction of emission still remains. 

A halo with the same kinematics as the disk (centre left 
panel) poorly reproduces the data: in this model the inner halo 
regions (/ e [0°, 90°] at vlos > - the 'receding' quadrant QI, 
and / G [-90°, 0°] at vlos < - the 'approaching' quadrant 
QIV) show 'horn-like' structures due to the fast rotation of the 
gas, whilst in the data these structures are absent. Also, the re- 
gions related to the external halo (/ g [0, 180°] at vlos < 
and / G [-180,0°] at vlos > 0) do not reach the same vlos 
of the data. Adding a vertical lag in the halo rotation (bottom 
left panel) improves the comparison with the data: the horns 
disappear thanks to the lower rotation, while the external halo 
reaches higher velocities along the line of sight. This compar- 
ison suggests that the H i halo of the Galaxy does not corotate 
with the disk. 

A vertical motion (central panels) globally shifts the halo 
emission to negative vlos for < or to positive vlos for > 0, 
because the gas, for any fixed b, is approaching the line of sight 
with an additional velocity sin(/?). This motion is not visible 
at low latitudes, but dominates at the highest ones. The high 
latitude emission in the LAB data occurs mostly at vlos < 
suggesting the presence of a global vertical infall of the gas, 
as already pointed out by the analysis of the HVCs and IVCs 
kinematics (Bajaja et al.|1987| ). A global vertical outflow from 
the disk is excluded. 

The efl'ects of the radial motions (right panels) are not eas- 
ily foreseen, because of the ofl'-centred position of the Sun with 
respect to the rotation axis. Obviously around the Anti-Centre 
region (/ = 180°) the emission is shifted toward negative (for 
Vr < 0) or positive (for vr > 0) velocities, because the gas is al- 
most directly approaching the Sun or receding from it. Around 
the Galactic Centre the l-v profile is globally 'deformed' (with 
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Fig. 2. l-v slices atb = 30° for the LAB dataset (top panel) and for six different halo models. Left panels: halo corotating with 
the disk (top) and rotating with vertical lag dV(p/d\z\ = -20kms"^ kpc"^ (bottom); central panels: halo corotating with vertical 
motion = -50kms~^ (top) and = +50kms~^ (bottom); right panels: halo corotating with radial motion vr - -50kms~^ 
(top) and Vr - +50kms"^ (bottom). Regions with |vdevI < 35 kms"^ have been excluded. All the cubes have been smoothed to 
8° resolution. Contour levels are: 40, 80, 160, 320, 640, 1280 mK in brightness temperature. 
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Fig. 3. The effect of the gradient- thickness degeneracy: the /-v slices at Z? = 45° for the two different rotating halo models are 
almost indistinguishable. The disk emission has been excluded and the cubes have been smoothed to 8° resolution. Contour levels 
are: 40, 80, 160, 320, 640, 1280 mK in brightness temperature. 



respect to the corotating case) toward positive velocities for 
V/? > and toward negative ones for vr < 0. 



This qualitative analysis suggests that the Hi halo of the 
Milky Way has two kinematical properties: lagging rotation 
and global infall. In principle we could quantify the various 



A. Marasco & R Fraternali: Modelling the Hi halo of the Milky Way 



5 



velocity components by comparing our models with the data. 
However these models have been built by assuming a value for 
the scale-height of the H i halo and we found this to be a crit- 
ical assumption as there is a degeneracy between the vertical 
rotational gradient and the halo scale-height. Because of pro- 
jection effects, halving hg and doubling dV(p/d\z\ (or vice versa) 
distributes the emission over the channels in similar ways at 
the high latitudes. As an example, Fig. [3] shows the l-v slice 
at = 45° for two halo models with circular rotation: the first 
with hg = 2kpc and dV(p/d\z\ = -20kms"^ kpc"\ the second 
with hg = 0.5 kpc and dV(p/d\z\ = -80kms"^ kpc"^ The total 
fluxes of the models have been normalized to the same value. It 
is evident that, even if the parameters used diff'er significantly, 
the first model is almost indistinguishable from the second. 

Since at high latitudes we are not able to discriminate be- 
tween a strongly lagging thin halo and a weakly lagging thick 
halo, we need an independent measurement for either hg or 
dV(p/d\z\. We tried to estimate hg from the emission along b in 
the H I velocity-integrated maps (which do not depend on the 
kinematics) but we found that this method is also degenerate. 
Consequently, we searched for an independent measurement of 
the velocity gradient. 

3. The rotation velocity gradient in the inner 
Galaxy 

3.1. Method of tangent-point 

The classical method to estimate the mid-plane rotation curve 
in the inner Galaxy (inside the solar circle) is that of the 
tangent-point: if the gas is in circular rotation without turbulent 
motions, for any direction / the highest velocity measured in the 
profile (the 'terminal velocity' Vter) is the line-of- sight velocity 
of the ring at the tangent-point radius /?tan = ^ol sin(/)| (the 'tan- 
gent velocity' Vtan). Therefore, using eq. ([T]) for /? = 0° we have 
V0(^tan, 0) = |vter + v© sin(/)|. The presence of turbulent motions 
smears the Hi line profiles, so Vter is actually lower than the 
highest velocity and it has to be searched in the wings of the 
profiles with some assumptions. Some authors used a bright- 
ness temperature threshold ( |Malhotra ]T995 ), others chose val- 
ues related to first peak in the spectrum ( Kerr 1962, 1964 j) or 
to the integral of the profile past this peak ( Shane & Bi eger-| 
|Smith|p^66| ). [Levine et al.| ( |2QQ8| ) used a more complex ap- 
proach based on the fitting of the whole line profile around the 
tangent velocity ( [Celnik et al.|1979 ). 

All those methods which make use of the end-line peak in 
the spectra - including that of Levine et al.| ( |2QQ8 ) - are not 
applicable at high latitudes since this peak, due to the falloff' 
in the Hi density at high z, is too faint to be detected. Thus, 
we adopted a diff'erent approach, often used for the external 
galaxies: we assumed a certain velocity dispersion for the gas, 
and, for each line profile, we performed a Gaussian fit to the 
emission falloff' leaving the amplitude and the mean (vter) of the 
Gaussian as free parameters. This method, if applied at the mid- 
plane, gives similar results as most of the approaches described 
above, but it also has the advantage to be suitable at diff'erent 
latitudes (see below). In exchange it requires a higher signal- 
to-noise (S/N) ratio with respect to other approaches, since it 



uses for each spectrum only the last points at low brightness 
temperature. To improve the S/N we smoothed the datacubes: 



tests performed on the toy models (Section 3.2) showed that 
after smoothing to 3° -4° this approach works well for values 
of S/N typical of the LAB data. 

To estimate the inner velocity gradient of the Milky Way we 
applied the tangent-point method outside the mid-plane. For 
any (/, b) in the inner Galaxy, the tangent point is at /?tan = 
7?o|sin(/)| and Ztan = cos(/) tan(Z7). However, if the gas is 
not exactly in corotation at any height from the plane, Vter is 
no longer equal to Vtan- For example, in a lagging halo for any 
/? 9^ 0° the highest velocity emission comes from a region closer 
to us than the tangent point, and if we use the method described 
we would assign to V0(7?tan,^tan) a value valid for a lower z: 
basically we would overestimate v^. The question is whether 
this geometric bias forbids the derivation of V0(7?, z), or there is 
a window in the (/?, z)-space where this eff'ect is negligible and 
therefore: 



^^(^tan.-^tan) 



Vter 



cos(/?) 



-h Vq sin(/) 



(4) 



We built 3D models to check this (Section [3^ and found that it 
is feasible to estimate the rotation up to a few kiloparsecs from 
the mid-plane. 

3.2. Tests using the models 

We applied the tangent-point method to diff'erent models to ver- 
ify that: a) in the presence of moderate gradients, the geometri- 
cal bias produces negligible errors and Jv^ /dz can be correctly 
derived; b) diff'erent density distributions weakly aff'ect the es- 
timate of the gradient. We built two pseudo-datacubes with the 
same Hi halo mass (5 x 10^ M©) and the same kinematics: 
v^(R,0) = 220 km s"^ and Jv^/JIzl = -20kms"^ kpc"^ In the 
first (model A) the Hi density distribution follows eq. ([2]), in 
the second (model B) the halo density decreases exponentially 
as \z\ increases and it is constant in R (see insets in Fig. [4]). 
The vertical scale-height is 2 kpc for both models. We added a 
velocity dispersion of 20 km s"^ for the halo gas. 



We applied the method described in Section 3.1 to both 
models, calculating V(p(R,z) for \z\ < 6 kpc using the input Hi 
velocity dispersion as the sigma of the Gaussian for fitting the 
profiles. In the absence of noise we are able to perform fits us- 
ing components with an arbitrary low Tb, taking into account 
only the dispersion falloff' of the lines. In this ideal case the 
resulting velocity fields of both models are perfectly recon- 
structed for 1 < 7? < 5 kpc and 1 < z < 4 kpc. Outside this window 



the geometric bias (see |3.1| ) is no longer negligible. 

We added a Gaussian noise to our models, with a dispersion 
of 50 mK (IcTrms of the Hanning- smoothed release of the LAB 
dataset), and we performed the fits using only the components 
with Tb > 2.5(Trms = 125 mK: the velocity field resulted unsat- 
isfactory for model B, showing that a higher S/N was required. 
To improve the S/N we smoothed the model cubes to 4° and 
performed a further Hanning Smoothing: this brought the aver- 
age IcTrms to 5 mK and the threshold of the fit to 12.5 mK. We 
re-performed the fits averaging all the profiles within the new 
beam and using only independent points, finding that for any 
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Fig. 4. Rotation velocity v^piz) (black points) averaged in the re- 
gion 3.5 < 7? < 5.5 kpc computed with the tangent point method 
for the model A (top panel) and B (bottom panel). The red thick 
line is the linear fit. The insets show the radial density distribu- 
tion from R = to Rq. 



fixed R between 2 and 5 kpc the V(p(z) curve could be fitted with 
a straight line in good agreement with the input slope. We av- 
eraged V(p(z) over a range of R in order to improve the statistics 
on the single measurements, at 4° of resolution there are about 
3 independent points for 3.5 < 7? < 5.5 kpc. 

Figure [4] shows the 7?-averaged V(p(z) curve derived for the 
two models with the method described. We also included a 
point for the mid-plane. We assumed that the error on the single 
measurement V(p(z) is given by the squared sum of three terms: 
a) the rms of the V0(z) distribution over 3.5 < R < 5.5 kpc; b) 
the maximum statistical error on the Gaussian fits computed in 
this distribution; c) a 'base' error of 4 km s"^ that takes into ac- 
count uncertainties in the gas velocity dispersion. Fitting the 
points with a straight line we obtained dV(p/d\z\ = -18.9 ± 
0.9kms~^ kpc~^ for model A and -21.6+ 1.9 km s~^ kpc~^ for 
model B, in good agreement with the input gradient. We tested 
further this method varying the vertical scale-height of model 
A from 0.75 kpc to 3 kpc, obtaining similar results. Our method 
depends weakly on the density distribution, so we can proceed 
to apply it to the LAB data. 

3.3. Application to the data 

The Hi distribution in the Milky Way is neither exactly ax- 
isymmetric nor /^-symmetric. Therefore we derived four V0(z) 
curves in the inner Galaxy (R<Rq), one for each quadrant (QI 
and QIV, see Section [23] ) and for the positive and the negative 
latitudes separately. We used the LAB dataset smoothed both 
spatially (4°) and in velocity (Manning), averaging the profiles 
within the new beam and using only the independent points. 
The threshold used is 2.5 x 6 mK for QI and 2.5 x 4 mK for QIV, 
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Fig. 5. Vertical rotational gradient in the inner Milky Way. Thin 
lines: rotation velocity V(p(z) averaged in the region 3.5 < R < 
5.5 kpc derived with the LAB dataset using the tangent point 
method. Each curve stands for a diff'erent region. Error bars are 
not shown. Black points: average V(p(z) curve. Thick red line: 
linear fit to the average curve. The yellow region stands for the 
2(T error on the gradient (Icr = 4kms"^ kpc"^). 



since the noises in regions at positive and negative declinations 
slightly diff'er. We assumed a velocity dispersion of 20kms~^ 
for the halo gas. This value has been estimated for NGC 891 



by |Qosterloo et al.| (|2QQ7 ). In the Milky Way, |Ford et aLl ( |2Q08 



found a cloud-to-cloud velocity dispersion of 18 km s~ for the 
Hi emission in the lower halo regions. The errors on the single 
fits were derived as described in Section IT2l 

We found that in the LAB Survey the statistical mean of the 
noise is not zero but instead systematically around -hlO mK. In 
the full resolution dataset the 2.50-^^^ threshold is well above 
these values, but after a 4° smoothing has been applied it be- 
comes of the same order. The nature of this emission is not 
clear, and systematic errors at this low level in the LAB Survey 
may exist (P. Kalberla, private communication). We preferred 
to re-set the noise zero point channel by channel, subtracting 
the average value estimated in regions not contaminated by 
the Galactic emission. This is equivalent to using for our fit 
a threshold higher than 2.5cr diff'erent for each channel map. 

Figure [5] shows the four 7?-averaged V0(z) curves derived 
for each quadrant, as well as the mean curve and its linear 
fit. We added two points belonging to the mid-plane, one for 
QI and one for QIV, averaged in 3.5 < 7? < 5.5 kpc. We esti- 
mated the rotation curve in the mid-plane V0(/?, 0) (see thick 
lines in Fig. [8]) using the LAB data smoothed to 2°, a thresh- 
old of 2 K and 2 km s~^ for the 'base' error. The linear fit gives 
-15 km s"^ kpc"^ for the average vertical falloff' in the rotation 
velocity. We estimated an error of 4kms"^ kpc"^ obtained us- 
ing as a constrain the point around z = 2 kpc; the shaded area 
in Fig. |5] shows the Icr confidence level. 

If vertical and radial motions are present they may aff'ect 
the calculation of V0(z), since we showed (Section |2.3[ ) that 
non-circular motions modify the shape of the l-v distribution. 
We checked this by adding a vertical infall (v^ = -40kms~^) 
and a radial inflow (vr = -30kms"^) to our models and re- 
estimating V(p(z). We found that in the presence of these motions 
the gradient would be overestimated in the QI regions and un- 
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derestimated in the QIV regions. However, the overall effect is 
inside the error bar of our gradient estimate, and we can neglect 
it. 

Recently it has been pointed out that the velocity of the 
solar circle v© could be signi ficantly h igher than the standard 
value of 220km s"^ dReid et al.|2009t . We note that the value 
of v© has no impact on our gradient estimate. Eq.|4] states that 
changing v© would shift up or down the rotation velocity at a 
certain longitude (i.e. at a certain R) for any z. For instance, 
in Fig. |5] all the points would be vertically shifted by the same 
amount, with no consequences for the fitted slope. 

Finally, we point out that the implicit assumption behind 
the method used is that in the regions investigated the gas 
is dense enough to show a significant emission. Testing the 
method with very diff'erent density distributions and scale- 



e) the external galaxies. We also disregarded components with 
I^LOsI > 250 km s"^ , since at these velocities only some isolated 
'Very High Velocity Clouds' contribute to the emission, and 
obviously our models do not attempt to reproduce these fea- 
tures. Furthermore, we excluded the QI and QIV regions for 
\b\ < 50° (the 'horn-like' structure, see Section [23] ) since their 
emission - for a fixed density distribution - gives information 
mainly on dV(p/d\zl which has already been set. 

4.2. The best model 

We found that the kinematical parameters that minimize 
the diff'erences between the models and the data are = 
-20kms"\ vr = -30kms"^ For the vertical scale-height we 
obtained 1.75 kpc computing the residuals using the squared 
sum and 1.5 kpc using the absolute sum, thus we set hg = 
1.6 kpc. The Hi halo mass above 2 disk scale-heights (see 
is 3.2 X 10^ M©. Table [l] summarizes the param- 



2.1 



heights (Section [3^ gave comforting results. However, if an 
unusually deep depletion of Hi were present in the inner 
Galactic regions our gradient would be overestimated. 

4. Results 

In the previous Section we derived the inner rotational gradient 
of the Milky Way Hi halo. We proceed now to estimate the 
remaining three free parameters of the halo: the vertical scale 
height (hg), the vertical component (v^) of the global motion 
and the radial one (vr). Given the inapplicability of the tangent- 
point method we can not derive the rotation velocity for R>Rq. 
Thus we assume that: a) the mid-plane rotation curve remains 
flat beyond Rq (see fig. [8|; b) the rotational gradient, derived 
for 3.5 < R < 5.5 kpc, is the same (-15kms"^ kpc"^) at each 
radius. 

4.1. Minimization of residuals 

We built a series of models with difl'erent combinations of the 
three parameters (hg, v^, vr) and searched for the values that 
minimize the difl'erences with the data. We explored the follow- 
ing ranges: from -lOOkms"^ to Okms"^ in steps of lOkms"^ 
for Vz, from -50kms~^ to 50kms~^ in steps of lOkms"^ for 
Vr, from 0.5 kpc to 4 kpc in steps of 0.25 kpc for hg(R). We 
excluded a positive range for since the data clearly do not 
show a systematic shift of the emission towards the positive 
velocities. For each model the total flux computed in the halo 
region was normalized to the observed halo flux by multi- 
plying each rB(/, Z?, vlos) for an appropriate factor, and since 
Tb oc A^Hi ^ M^^^ this is a correction to the Hi halo mass. The 
resulting halo masses above 2 disk scale-heights for difl'erent 
models vary from 0.5 x lO^M© to 9 x lO^M©. 

For each model we computed the difl'erence with the LAB 
datacube pixel by pixel, and we summed these values using 
a cos(Z7) factor to take into account the projection efl'ects. We 
used both a quadratic sum and an absolute sum, obtaining lit- 
tle difl'erences. The sum were extended to all the components 
with the exception of the regions occupied by the disk emission 

(IvdevI < 35kms~^) and by: a) the Magellanic Stream; b) the the analysis described in Section \4A 



Section 

eters of the best halo model, as well as their errors. These latter 
were obtained by searching for 'acceptable models' around the 
best set of parameters. Note that the H i halo mass reported cor- 
responds to 5-10% of the total Galactic H i. 

Figure [6] shows seven representative / - v slices at various 
latitudes for three datacubes: a model with only the Hi disk 
component, a model with both (disk -h best halo) components, 
and the data. The H i distribution in the disk was modelled as 
reported in KD08, but we used the rotation curve found in 
Section 33 (thick lines in Fig. [8|, we added an isotropic ve- 
locity dispersion of lOkms"^ and we included an exponen- 
tial decrease in density for R < 4.75 kpc. We did not include 
the warp (but see Section [43] ). The comparison between mod- 
els and data clearly highlights how the H i disk emission alone 
fails to reproduce the observations, especially at the intermedi- 
ate latitudes (Z? = ±15° and b = ±30°). Instead, adding the halo 
component allows us to reproduce most of the H i features at 
low brightness temperature (Tb<0.5 K), especially those com- 
ing from the external regions (0° < / < 180° at vlos < 0. and 
-180° < / < 0° at Vlos > 0) at latitudes \b\ > 15°. The inner re- 
gions (QI and QIV) appear to reach velocities higher than the 
data: this can be interpreted as evidence of a larger gradient 
in the innermost regions - as found by [Qosterloo et al.| ( |2007 



for NGC 891 - or of a depletion of Hi in these areas (see next 
Section). Also, the LAB data shows asymmetries between the 
positive and the negative latitudes that we can not reproduce 
with our symmetric models. 

In order to estimate the impact of the local Hi emission 
on the above parameters we performed the following test. We 
considered three models, each one excluding a difl'erent part of 
the halo: in model #1 we included only the inner halo by ex- 
cluding the H I emission at 7? > 8 kpc, in model #2 we added 
the local gas by excluding the emission at 7? > 9 kpc, in model 
#3 we excluded the H i at 7? < 9 kpc in order to take into ac- 
count only the outer emission. For each model we performed 

We found that the best 



Leading Arm; c) the GCP complex (or 'Smith Cloud'; Smith 
|1963| ) which is an isolated cloud with an anomalous kinemat- 
ics ( [Lockman et al.|20Q8] ); d) the Outer Arm (see Section [2J] ); 



parameters for model #2 are the same as those of table [T] 
while those of model #1 and #3 slightly difl'er: we obtained 
(v^ = -25kms"^v/? = -30kms"\/Zg = 1 . 1 kpc) for the first 
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Fig. 6. /-V slices at 7 different latitudes indicated at the top left corners of the rightmost plots. First column', the Hi disk model; 
second column: Hi disk + best halo model; third column: the LAB data. Each datacube is smoothed to 8° resolution. Contour 
levels in brightness temperature range from 0.04 K to 81.92 K scaling by a factor 2. 
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Fig. 7. l-v slices at /? = +15° (top) and b = -l5° (bottom). First column: the warped H i disk model; second column: the warped best 
model (Hi disk + halo); third column: the LAB data. Each datacube is smoothed to 8° resolution. Contour levels in brightness 
temperature range from 0.04 K to 81.92 K scaling by a factor 2. 
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4.3. Further refinements 

We considered the effect of a Galactic warp. Following 'Levine 



Table 1. Global parameters of the H i halo of the Milky Way, 
from our best halo model (see text). 



and (v^ = -15kms~\v/? = -20kms"\/zg = 1.6kpc) for the 
latter. These values are all consistent within the errors with 
those of table[T]and we can conclude that our results are global. 

To understand if a variation of the parameters Rg and y af- 
fects our results, we build four new models using respectively 
7?g = 2.4 kpc, = 1.1 kpc, 7 = 2 and 7 = 8 (see eq.|2]). For 
each model we re-performed the residuals minimization, find- 
ing that in all cases the best parameters are consistent within 
the errors with those of table [T] Finally, we investigated the 
impact of v© on the residuals by re-building all models using 
Vq = 250 km s"^ for both the mid-plane rotation and the ve- 
locity projections, obtaining no variation from the above best 
parameters. 



et al.| p006 ), we modeled the off'set of the mid-plane zo from 
the plane z = with the formula 

Zo(R,(p) = Wo(7?)+Wi(7?)sin(0-0i(7?))+W2(^)sin(20-02(^))(5) 

where the dependencies on R of the three amplitudes Wi and 
the two phases cpj have been taken fromlKalberla et al.|(|2007| 



Fig. 16 and Fig. 17 therein). We considered the warp only for 
R > Rq, and we set W2(R) = Okpc for R < 16 kpc. Under 
the assumption that the H i halo stratifies like the H i disk, we 
included the warp in our models just substituting the density 
computed in (R, z) with that computed in (R, z - Zo) for both 
the disk and the halo distribution. Also, the halo rotation curve 
V(p(R, z) has been substituted with V(p(R, z - Zo). This can not be 
considered a full treatment of the warp effect since we used eq. 
^ to compute the line of sight velocity and we assumed that 
the orbits are circular. However, the mid-plane significantly de- 
viates from the z = plane only at large radii (R > 15 kpc) 
and around ~ 80°, which makes the contribution of the halo 
emission in the warp region negligible. 

Figure [T] shows the two / - v plots atb = +15° for the data 
(rightmost panels) compared with a disk-only model (left col- 
umn) and a disk+halo model, both with the inclusion of the 
warp. The warp is very asymmetric and mostly visible at pos- 
itive latitudes for / = 90° . In this region - the Outer Arm zone 
- the H I emission is deformed towards negative velocities and 
the model reproduces the data significantly better. At higher 
latitudes the effect of the warp is negligible and the other / - v 
diagrams in Fig.[6]remain unchanged. Clearly, the inclusion of 
a warp does not significantly affect our results for the H i halo. 

In order to reproduce the emission related to the inner re- 
gions (QI and QIV) we considered two effects separately. The 
first is the inclusion of a radial dependence for the rotational 
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Fig. 9. l-v slices at Z? = +15° (top) and = -15° (bottom) at longitudes -90° < / < 90°. First column: the best model shown 
in Fig.|6j second column: the model refined including a radial dependence of dV(p/d\z\ for R < 4.5 kpc; third column: the model 
refined including a linear decrease of the halo thickness for R<Rq\ fourth column: the LAB data. Each datacube is smoothed to 
8° resolution. Contour levels in brightness temperature from 0.04 K to 81.92 K scaling with a factor 2. 
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Fig. 8. Rotation curves as a function of the height from the 
Milky Way mid-plane used in modelling. The thick lines on 
top stand for the mid-plane rotation. Red solid lines: reced- 
ing curves; blue dashed lines: approaching curves; black dotted 
lines: rotation velocities assumed for the 'shallow rise' model 
at 7? < 4.5 kpc. For R>1 .1 kpc we assumed that the curves are 
flat. 



gradient in the inner halo. We assumed that for R < 4.5 kpc 
the magnitude of the rotational lag increases linearly each kilo- 
parsec by lOkms"^ kpc~^ The efl'ect of this 'shallow rising' 
rotation curve (dotted curves in Fig. [5]) is to lower vlqs in the 
inner halo, reducing the 'horns' (see Section [23] ) in the l-v dia- 
grams. Figure |9] shows that this eff'ect is not sufficient, since the 
l-v profiles obtained for this model (second column) diff'er only 
slightly from the one with a constant vertical velocity gradient 
(first column). 



As a second attempt we reduced the halo thickness in in- 
ner regions. We assumed that the vertical scale-height of the 
Hi halo decreases linearly for R<Rq, vanishing at the Galactic 
Centre. Figure [9] shows a model with this thin inner halo (third 
column) and constant rotational gradient of -15 km s"^ kpc"^ . 
This model reproduces fairly well the H i emission in these re- 
gions of the datacube. 

We conclude that the scale-height of the Hi halo of the 
Milky Way inside Rq is likely to be lower that our average 
value of 1.6 kpc, the rotational gradient may also be larger than 
15 kms"^ kpc"^ in the inner parts. It is difficult to disentangle 
these two eff'ects in our analysis. For instance, given the de- 
generacy between dV(p/dz and hg (Section [23] ), a non-constant 
gradient with R would imply a corresponding variation of the 
vertical scale-height. Both a decreasing scale-height and an in- 
creasing gradient in the inner parts are consistent with a galac- 
tic fountain origin for the Milky Way halo (see Section [5?T] ). 

5. Discussion 

In this paper we have estimated the global properties of the 
Milky Way Hi halo by building 3D models and comparing 
them with the LAB Survey dataset. We found that the Galactic 
Hi halo can be globally described by the following average 
kinematical parameters: dV(p/d\z\ = -15 ± 4k ms~^ kpc~\ = 
-20+^kms"^ and vr = -30+^kms"^ (Section 
lowing we discuss the results obtained. 



4.2). In the fol- 



5.1. Physical interpretation 



Fraternali & Binney[ ( |2006| ) successfully reproduced the verti- 
cal gaseous distribution of NGC 891 and NGC2403 by mod- 
eling the halo as a continuous flow of non-interacting galac- 
tic fountain clouds moving ballistically in the galactic poten- 
tial (see also [Collins et al.||2002t[Heald et al.|2QQ7] ). However 
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they were forced to include interactions of these clouds with 
the environment in order to reproduce the vertical rotational 
lag ( [Fraternali & Binney| |2008 ): in this scenario the fountain 
clouds accrete gas from an ambient medium with a low an- 
gular momentum with respect to the galactic disk. This process 
slows down the rotation of the galactic fountain clouds and pro- 
duces gas accretion towards the disk. Hydrodynamical simula- 
tions show that the Galactic corona could be the reservoir of 
this accreting gas f Marinacci et al.||2Q10a ). The typical orbits 
of this 'accreting galactic fountain' differ from those of a foun- 
tain without accretion for the presence of inward radial motions 
(Fig. 6 in Fraternali & Binney||20Q8 ) which indeed have been 
observed in the Hi halo of NGC 2403 ( [Frate maH et al.|2QQl|). 
In Section 4.2 we showed that the H i halo of the Milky Way has 
global vertical and radial motions of inflow, and in Section 3.3 
we estimated a vertical lag equal to the one found in NGC 891. 
These results suggest that the dynamical interpretation given 
for the gaseous halo of NGC 891 is suitable for the Milky Way 
halo as well. 

The vertical inflow motion estimated (v^ < 0) can be ex- 
plained supposing that the galactic fountain gas is initially ion- 
ized for a certain part of the rising orbit, becoming visible in 
H I only after recombination. If one considers a typical foun- 
tain orbit where the gas is ejected vertically from the disk with 
an ejection velocity Vkick and a vertical velocity that decreases 



roughly linearly with time (Fig. 6 in |Fratemali & Binney 2008 ), 
to have an average vertical velocity < the gas should be 
ionized for a fraction -2vJ/vkick of the rising orbit. This frac- 
tion is around 50 - 70% for the vertical velocity derived in 
^) and for typical ejection velocities 



4.2 



(-20+^kms" 



Section 
of 60-80kms"\ 

The negative value of vr in the halo is more difficult to in- 
terpret. Note that, in the context of a phase-change fountain, 
vr<0 does not imply a net inflow for the whole gaseous halo 
but only for the neutral descending clouds. For each cloud, this 
inward motion would be roughly balanced by an outward mo- 
tion during the first part of the orbit, when the gas is ionised. 
A ballistic non-interacting fountain cloud may have some nega- 



tive Vr, but only at the very end of its orbit ( [Fratemali & Binney 
2006| ). It seems therefore more likely that such a large value of 
Vr is produced by the loss of angular momentum due to gas ac- 



cretion hypothesized by Frate rnali & Binney | ( |2008| ). In order 
to fully understand these effect, one requires a full dynamical 
modelling of the halo, which will be the subject of further in- 
vestigations. 

5.2. Relation to previous works 



Lockman (1984]) found that in the lower halo regions of the 
inner Galaxy a large amount of Hi is present. He estimated 
that this medium had a vertical scale-height of 500 pc approx- 
imately constant for 4 < 7? < 8 kpc. Also, iLockman] ( |2002| ) 
confirmed that the diffuse H i emission in these areas is orga- 
nized into discrete clouds (see also |Ford et al.|2008] ). In Section 
4~2| we found that all the acceptable models have larger scale- 
heights. The difference can be explained as follows: a) our av- 
erage scale-height is derived using the emission from the entire 



halo while it seems in the innermost regions of the Galaxy the 
halo is thinner (Section [4~2| and [43] ); b) Lockman derived the 
position (R, z) of the emission at the tangent point assuming 
co-rotation for the halo gas, and he pointed out that in case of 
a lagging halo the amount of H i would be underestimated. 



[Levine et al.[(|2008]) were the first to derive a value for the 
vertical falloff in the Galactic rotation velocity using a direct 
approach. They applied the tangent point method to both the 



VLA Galactic Plane Survey (StiletaL 2006) and Southern 



Galactic Plane Survey ( McClure-Griffiths et a 



] [2QQ5] ) dataset 



to derive the inner rotation curves up to |z| = 100 pc. The method 
used to determine the tangent velocity was based on fitting the 
whole profile around Vter using a form derived by [Celnik et al.[ 
(1979). They found an average lag of -22 ± 6kms~ kpc" , 
in agreement with our estimate obtained in the halo region for 
1 < |z| < 4 kpc. With the resolution and the S/N of the LAB data 
we can not determine the rotation curve for |z| < 1 kpc, thus a 
direct comparison with the work of Levine et al. is not possible. 

A local value for the vertical rotational gradient was derived 
by Pidopryhora et al. [ ([2007[) in a small region (the 'plume') 
above the Ophiucus superbubble. Assuming that the feature is 
located at the tangent point, they found a lag of -27 km s"^ with 
respect to the corotation at position R = 4 kpc, z = 3.4 kpc. This 
means an average gradient of -8 km s"^ kpc"^ at that radius. In 
the same region we measured a gradient of roughly a factor 2 
larger (see receding northern curve in Fig.[5|. In the LAB data 
smoothed at 4°, the Ophiucus superbubble is not visible and we 
are likely measuring the average velocity of the surrounding 
region where the superbubble is expanding. 

In the attempt to reconstruct the spatial H i distribution by 
a 3D de-projection of the LAB data, rotation curves at several 
heights from the mid-plane have been derived by [Kalberla et| 
^aD (2 001\ usin g a theoretical approach, previously adopted by 
Levine et al.[ ( [2006[ ). They modelled the Galactic dark mat- 
ter distribution so that the resulting rotation curve (obtained 
assuming that the gravity dominates the gas dynamics every- 
where) gives an hydrostatic Hi distribution whose properties 
(e.g. the flaring) are consistent with the dark halo itself. The 
curves obtained (Fig. 9 in [Kalberla et al.[2007j ) differ from those 
we derived (Fig. [8]), specifically the vertical lag at = 4.5 av- 
eraged for 1 < |z| < 4 kpc is around -22kms~^ kpc~^ and it 
decreases as R increases. Our approach does not rely on any 
assumption about the dynamical state of the system as we sim- 
ply derive the best model that reproduces the observations. We 
notice that if the H i halo is produced by galactic fountains in- 
teracting with the environment (see Section [5^] ), assuming that 
the dynamics of the gas is dominated by gravity is not neces- 



sarily a good description of the system (Fratemali & Binney 
[2008] ). [Maririacci et al.[ ( [20TQb] ) also showed that for NGC 891 
a hydrostatic model does not reproduce satisfactorily the data. 



model derived in 
10% of the Hi is not 
related to the disk, in agreement with our estimate for the H i 
halo mass (Section [42]). 



FinaUy, KD08 analysed the Galactic H 
Kalberla et aT] ( [2QQ7] ). They found that 



12 



A. Marasco & F. Fraternali: Modelling the Hi halo of the Milky Way 



5.3. HVCs and IVCs 

HVCs and IVCs are regarded as the main evidence for the pres- 
ence of a gaseous halo in our Galaxy, and it is important to un- 
derstand how they fit with the results of this work. Even if the 
emission of these objects dominates the H i flux in regions out- 
side the Hi disk, at low column densities the sky is basically 
filled with neutral hydrogen emission (Fig.[T]). 

Distance estimates for the IVCs show that they are a local 
phenomenon ( Wakker||2Q04 and references therein) and their 
disk-like metallicities and kinematics suggest a galactic foun- 
tain origin. In this paper we showed that these clouds are em- 
bedded in a diff'use halo, which shares the same kinematics. 
The emission of this gaseous component (IVCs -h diff'use H i) 
is reproduced by our model. This suggests that the whole H i 
halo is made up of IVC-like objects, most of which are at large 
distances and thus unresolved. If these clouds had an average 
H I mass of ~ lO^M©, the H i halo would contain 3-4 thousand 
such objects. 

Recent studies on HVCs established that these objects have 
higher distances and heights from the plane and lower metal- 
licities than IVCs ( jWakker et aLl|2QQ7| |2QQ8| ). In general our 
halo model does not reproduce HVCs and the halo parameters 



derived in Section [42] are not representative of their properties. 
This suggests that the nature and origin of these clouds are dif- 
ferent from that of the halo described here. There are special 
cases such as the the well-known complexes C and A, visible 
in Fig.[6]at b = -h50° and vlos ~ -140 km s"^ and b = -h30° and 
vlos — 160 km s"^ Although separated from the 'normal' halo 
gas these clouds appear in a location that is not completely pe- 
culiar. They follow the trend of the halo emission but at a more 
extreme velocity. A recent study puts complex C at ~ 7 - 8 kpc 
above the mid-plane ( Thom et al._2008 ). Interestingly, its pro- 
jected velocity is not very different from what one would de- 
rive by extrapolating our vertical rotational gradient at these 
heights. 

5.4. Comparison with external galaxies 

The properties of the Galactic H i halo, derived here, are simi- 
lar to those of external galaxies. The vertical rotational gradient 



estima ted in Section 3.3 is the same found by Oosterloo et al. 
( |2007t for NGC 891; inward radial motions of 10 - 20 km s"' 
have been found by [Fraternali et al. (2002) in the Hi halo of 

poos] )' showed the presence of ver- 



NGC 2403; Boomsma et al. 



tical motions for halo gas clouds in NGC 6946. The halo Hi 
mass derived for the Milky Way is also consistent with the val- 
ues estimated in other galaxies (see [Fraternali|20 1 0) . However, 
we notice that the techniques used to separate halo gas from 
disk gas are very diff'erent and a direct comparison is difficult. 
Nevertheless, it is clear that a remarkably massive halo like that 
detected in NGC 891 (1.2 x lO^M©, 30% of the total Hi) is not 
present in the Milky Way and it is perhaps fair to say at this 
point that NGC 891 constitutes an extreme case. 

A recurring question is whether or not the anomalous H i 
clouds observed around external galaxies are of the same na- 
ture as the galactic IVCs and HVCs. If all halos are made up 
of clouds with properties similar to the IVCs then in external 



galaxies they would be unresolved and appear like a smooth 
medium, as is observed. Observations of external galaxies 
also show the presence of very massive clouds (~ lO^M©) 
such as the filaments of NGC 891 ( [Oosterloo et al. [12007 ) and 
NGC 2403 ( [Fraternali erar||2Q02l ). These clouds have similar 
properties to large Galactic HVCs like complex C. They have 
the same masses (the Hi mass of complex C is ~ 5 x 10^ M©, 
[Thom et al.|2008| ) and they show comparable kinematics. For 
instance the filament in NGC 2403 follows the general kine- 
matic pattern of the halo gas but at a more extreme velocity. 
This is analogous to what has been described in Section 53_ 
for complex C. In external galaxies, clouds are also detected 
at large distances and/or very anomalous velocities: e.g. the 
'forbidden' clouds in NGC 2403 (Fratern ali et al.[|2QQ2l ), the 
so-called counter-rotating clouds (see Fratemali|2010| ) and the 
clouds surrounding M 3 1 at tens of kpc from the disk ( [Thilker 
et al.|2004| ). Very massive and very anomalous clouds are not 



consistent with a galactic fountain origin. They are more likely 
to be gas falling into these galaxies for the first time either from 
satellites or directly from the IGM, although the latter mecha- 
nism is not clearly understood ( [Binney et al.|2009| ). All things 
considered, they are similar to the Galactic HVCs. 



6. Conclusions 

Our Galaxy off'ers the unique opportunity to study an H i halo 
from two diff'erent perspectives: on the one hand one may anal- 
yse the single bright anomalous clouds - the HVCs and IVCs - 
which populate the extra-planar regions, on the other hand one 
may consider the halo as a global component and analyse its 
global parameters. In this paper we adopted the latter approach. 

We modelled the extra-planar Hi of the Milky Way 
and compared its emission with the Leiden- Argentine-Bonn 
Survey. Using the tangent point method we estimated the ver- 
tical gradient in the rotation velocity. The main results can be 
summarized as follows: 

1 . the Milky Way has a halo of neutral gas with an H i mass 
of 3.2+^-^ X 10^ Mo (~ 5-10% of the total Hi mass) and an 
average vertical scale-height of 1.6^^4 kpc, which seems to 
be decreasing in the innermost regions. 

2. The Hi halo has a global vertical lag in rotation velocity. 
We measured the magnitude of this lag in the region 3.5 < 
R<5.5 kpc and 1 < |z| < 4 kpc, obtaining an average value of 
-15 ± 4kms"^ kpc"^ 

3. We detected a global inflow motion in both the vertical 
(-20^^ km s"^) and the radial (-30^^ km s"^) directions. 

4. The Hi halo is likely to have been generated by super- 
nova explosions, which create a galactic fountain. If a sin- 
gle fountain cloud has properties similar to the local IVCs, 
then the halo should contain thousands of these objects. 

5. The properties of the Milky Way halo are broadly similar 
to those derived for external galaxies. 
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